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Abstract (Basic) : EP 624847 A 

The generator produces an approximation function based on pairs of 
values (Xi,Yi) associating a dependent magnitude (Yi) with an 
independent magnitude (Xi) , and generates new data pairs (Xa,Ya) 
following the approximation. 

The generator has a first part (10) that iteratively determines a 
current linear regression function to select those of the current 
regression functions that deliver the approximation to the current set 
of input values with least error. The parameters (p,q) of the best 
regression function are then retained, and used by the second part (17) 
to generate the new approximate pair of data values. 

ADVANTAGE - Requires little computational capacity and gives rapid 
results under controlled maximum error, for use in neural networks or 
in approximation to measured values in industrial processes. 
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Abstract (Equivalent) : US 5519647 A 

An apparatus (5) for generating an approximation function based on 
first pairs ((XI, Yl) to (X6, Y6) ) of values associating a dependent 
variable (Yl to Y6) with an independent variable (XI to X6) , and for 
determining second pairs (XA, Y ' A) of values of said variables in 
accordance with said approximation function, characterized in that the 
apparatus comprises: 

first means (10) : 

for iteratively determining at least one current linear regression 
function by making first errors of alternating sign equal in absolute 
value, which first errors have been measured between, respectively, 
first values (Y3, Y4, Y5) of the dependent variable for three pairs 
(X3, Y3) (X4, Y4) (X5, Y5) of a series of pairs, and second values 
(Y f 3, Y'4, Y*5) of the dependent variable determined, in accordance 
with said current linear regression function, for the same values (X3, 
X4, X5) of the independent variable, 

for selecting that one of the current linear regression functions 
which produces an approximation of all the pairs of said series with a 
minimum error, and 

for coding the selected linear regression function with the aid of 
specific codes (p, q) , and 

second means (17) for determining said second pairs (XA, Y 1 A) with 
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© Dispositif (5) pour generer une fonction d'approximation fondee sur des premiers couples ((Xi. Yi) a (Xs. 
Y 6 )) de valeurs associant une grandeur dependante (Yi a Y 6 ) a une grandeur independante (X, a Xg). et pour 
determiner des seconds couples (X A( Y' A ) de valeurs desdites grandeurs d'apres ladite fonction d'approximation. 
Le dispositif comprend : 

- des premiers moyens (10) pour determiner iterativemant au moins une fonction lineaire courante de 
regression, pour selectionner celie des fonctions lineaires courantes qui delivre I'approximation de tous les 
couples de ladite suite avec des erreurs minimales, et pour coder, a I'aide de codes specifiques (p, q), la 
fonction lineaire de regression selectionnee. 

- et des seconds moyens (17) pour determiner lesdits seconds couples <X A , Y* A ) a I'aide desdits cod s 
specifiques. 

Le dispositif peut etre utilise pour le calcul de valeurs approchees de fonctions mathematiques, par exemple 
dans un reseau de neurones, ou pour determiner une fonction de regression formant une approximation de 
valeurs experimentales de mesure. par exemple des mesures dispersees issues du suivi d'un processus 
industriel. 

L'invention concerne egalement une methode pour generer une fonction d'approximation. 
Application : calcul approche de .fonctions mathematiques, calcul de fonctions de regression. 
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L'invention concerne un dispositif et une methode pour generer une fonction d'approximation, celle-ci 
etant fondee sur des premiers couples de valeurs associant une grandeur dependante a une grandeur 
independante, et pour determiner des seconds couples de valeurs desdites grandeurs a partir de tadite 
fonction d'approximation. 

5 Un dispositif et une methode de ce type sont connus du brevet US-A- 3 789 203 qui decrit un 
generateur de fonctions operant une approximation par interpolation iterative. Ce dispositif est prevu pour 
etre utilise dans des applications de traitement de donnees necessitant un calcul de fonctions telles que 
sin(x). tg(x) par exemple. Ce dispositif ne requiert qu'une capacite de stockage minimale de la part d'un 
dispositif utilisateur. A partir de deux points appartenant a une fonction a interpoler. la methode tout d'abord 

w interpole la fonction par une droite reliant les deux points, puis fait une approximation des ecarts entre la 
droite et la fonction par des approximations polynomials de degre croissant. Ensuite eile substitue aux 
points initiaux des points approximatifs pour reduire la longueur du segment reliant les points a traiter et 
enfin reitere les operations precedentes. 

Une telle methode necessite des ressources importantes en moyens de calcul et ne peut etre mise en 

is oeuvre qu'avec des calculateurs performants. 

Or, il existe des applications ou une teHe methode ne peut pas etre mise en oeuvre car elles ne 
disposent pas des ressources suffisantes. De plus, pour certaines applications on peut se satisfaire d'un 
calcul approche de la fonction pour des valeurs, en nombre limite. de la grandeur independante. 

II peut s'agir d'une fonction sigmoVde appliquee a des potentiels neuronaux delivres par au moins un 

20 neurone dans un reseau de neurones. II peut s'agir d'une autre fonction non lineaire, par exemple une 
fonction racine, pour des calculs de distances entre des etats de neurones. Les applications peuvent aussi 
concerner d'autres dispositif s com me des generateurs de fonctions, des calculateurs ou autres. 

Pour calculer une telle fonction, sans passer par une fonction d'approximation, on peut utiliser 
differentes manieres. 

25 On peut effectuer le calcul mathematique exact pour chaque valeur de la grandeur independante a 
traiter, en programmant un calculateur selon les methodes connues. Une telle methode necessite d'effec- 
tuer a chaque fois les memes operations ce qui peut necessiter beaucoup de temps si le nombre de 
valeurs est eieve. 

On peut aussi prealablement stocker dans une memoire des tables precalculees. Dans ce cas. la 
30 lecture en memoire du resultat peut etre rapide. Mais pour couvrir avec un pas assez fin toutes les valeurs 
possibles de la grandeur independante, il faut alors disposer de tables de grandes capacites. Ces methodes 
de calcul presentent done des inconvenients. 

D'autre part, on peut etre conduit a identifier deux grandeurs qui sont dans la dependance Tune de 
I'autre par des couples de valeurs associant une grandeur dependante a une grandeur independante. Ainsi, 
35 dans le suivi d'un processus industries on peut etre conduit a mesurer par exemple un rendement R d'une 
operation en fonction de la temperature T a laquelle a ete realisee ladite opeation R = f(T). Un suivi du 
processus peut consister a consigner des lots de couples de mesures sur un graphique. Celui-ci peut etre 
utilise pour car ac tenser le processus ou pour en deduire de nouveaux parametres de regiage de ladite 
operation. Ceci fait par exemple I'objet de I'article de H. ISHIBUCHI et H. TANAKA, "Regression analysis 
40 with interval model by neural networks" in "IEEE International Joint Conference on Neural Networks", vol. 2, 
18-21 Nov. 1991, SINGAPORE. II est souhaitable que ces nouveaux parametres soient representatifs du 
fondement de ladite operation et que les fluctuations de mesures inherentes a ce genre de processus 
soient ecartees. II est done souhaitable de determiner une approximation de ta fonction f{.). 

Ainsi dans un cas il peut s'agir de mesures erratiques ou entachees d'erreurs que Ton desire 
45 representer par une fonction d'approximation. 

Dans un autre cas, on connait des valeurs precises mais I'utilisation a en faire ne necessite pas une 
grande precision et une fonction d'approximation sufftt. 

Un des buts de ('invention est de generer une fonction d'approximation avec des moyens materiels 
reduits permettant de calculer rapidement un nombre reduit de valeurs de la grandeur dependante utiles a 
so I'application sans avoir pour cela a determiner d'autres valeurs de la fonction d'approximation. Un but 
complementaire est de delivrer des valeurs qui peuvent etre approchees dans la limite d'une erreur 
maximate controlee. 

Ce but est atteint avec un dispositif caracterise en ce qu'il comprend : 
- des premiers moyens : 

55 - pour determiner iterativemant au moins une fonction lineaire courante de regression en rendant 

egales, en valeur absolue, des premieres erreurs de signes alternes mesurees entre des premieres 
valeurs de la grandeur dependante pour trois couples d'une suite desdits premiers couples, et 
respectivement des secondes valeurs de la grandeur dependante determinees, d'apres ladite 
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fonction lineaire courante. pour les memes valeurs de la grandeur independante, 

- pour selectionner celie des fonctions lineaires courantes qui delivre I'approximation de tous tes 
couples de ladite suite avec des erreurs minimales. 

- et pour coder, a I'aide de codes specifiques, la fonction lineaire de regression seiectionnee, 

5 - et des seconds moyens pour determiner lesdits seconds couples a I'aide desdits codes specifiques. 

Ainsi avantageusement on determine une fonction lineaire de regression approchant au mieux les 
differents couples de valeurs connus. Les resultats approches ainsi delivres forment un compromis 
sattsfaisant pour de nombreuses utilisations du dispositif generateur de fonctio... 

Une fonction lineaire de regression est une fonction simplificatrice qui represente un phenomene 
w complexe en reduisant les parametres significatifs. En representant la suite de couples de valeurs par des 
points dans un espace a deux dimensions, la fonction lineaire de regression devient une droite de 
regression. 

Ainsi apres avoir defini la droite de regression par des codes, on peut calculer une valeur approchee de 
la grandeur dependante en tout point de la droite de regression avec des moyens reduits pour des valeurs 
;5 quelconques de la grandeur independante. 

L'invention concerne egalement une methode pour generer une fonction d'approximation. la methode 
comprenant : 

- une premiere phase : 

- pour determiner iterativement au moins une fonction lineaire courante de regression en rendant 
20 egales, en valeur absolue, des premieres erreurs. de signes alternes. mesurees entre des 

premieres valeurs de la grandeur dependante pour trois couples d'une suite desdits premiers 
couples, et respectivement des secondes valeurs de la grandeur dependante determinees, d'apres 
ladite fonction lineaire. pour les memes valeurs de la grandeur independante, 

- pour selectionner celle des fonctions lineaires courantes qui delivre ['approximation de tous les 
25 couples de ladite suite avec des erreurs minimales. 

* et pour coder la fonction lineaire de regression seiectionnee a I'aide de codes specifiques. 

- et une seconde phase pour determiner lesdits seconds couples a I'aide desdits codes specifiques. 
Les moyens mis en oeuvre par l'invention peuvent etre formes par un calculateur programme ou par un 

circuit dedie. lis peuvent aussi mettre en oeuvre des neurones. 

30 Un dispositif mettant en oeuvre des neurones selon l'invention peut etre utilise par un reseau de 
neurones, dont il peut notamment en constituer un sous-ensemble. En effet, pour fonctionner, le reseau de 
neurones doit disposer de moyens pour appliquer une fonction non lineaire d'activation aux potentiels de 
neurones qu'il delivre. Selon l'invention, le dispositif muni de neurones peut calculer une approximation de 
cette fonction non-lineaire d'activation. II peut egalement calculer des distances entre des etats de neurones 

35 en calculant une approximation d'une fonction racine carree destinee a etre exploitee dans le reseau de 
neurones. 

Lorsque la taille de la suite de couples de valeurs fournis initialement est elevee. on peut diviser la suite 
de couples en plusieurs sous-ensembles pour determiner plusieurs droites de regression et ameliorer la 
precision de 1'approximation. La fonction d'approximation de la suite de couples est alors formee par une 
40 fonction lineaire par morceaux pour laquelle une exigeance de continuite entre les morceaux peut etre ou 
non imposee. 

Certains couples de ('ensemble de couples de valeurs peuvent avoir une influence particuliere que Ton 
peut concretiser en donnant un coefficient de ponderation specifique a chaque couple. Dans ce cas. I'erreur 
affectee a chaque couple tient compte de ce coefficient de ponderation specifique. 
45 Ces differents aspects de l'invention et d'autres encore seront apparents et elucides a partir des modes 
de realisation decrits ci-apres. 

L'invention sera mieux comprise a I'aide des figures suivantes donnees a titre d'exemples non limitatifs 
qui representent : 

Figure 1 : un graphique montrant une representation a deux dimensions d'un ensemble de points avec 
so une droite de regression D. 

Figure 2 : un graphique montrant un ensemble de points et des droites servant a la determination d'une 
enveloppe. 

Figure 3 : un organigramme d'une premiere variante de mise en oeuvre de la methode a partir de 
triplets de points. 

55 Figure 4 : un organigramme d'une seconde variante de mise en oeuvre de la methode a partir de triplets 
de points. 

Figure 5 : une partie d'organigramme d'une troisieme variante de mise en oeuvre de la methode a partir 
de couples de points. 
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Figure 6 : une partie d' organigramme indiquant une preselection de points appartenant a une enveloppe 

inferieure ou a une enveloppe superieure de I'ensemble des points. 

Figure 7 : un organigramme pour la determination des enveloppes inferieure et superieure. 

Figure 8 : un schema d'un dispositif selon ['invention. 

Figure 9 : un schema d'un dispositif de codage pour le calcul et le codage d'une droite de regression 
dans le cas general. 

Figure 10 : un schema d'un dispositif de codage pour le calcul et le codage d'une droite de regression 
pour la premiere variante. 

Figure 1 1 : un schema d'un dispositif de transcodage pour le calcul des valeurs de la grandeur 
dependante a partir d'un ensemble de droites de regression codees. 

Figure 12 : un schema d'une realisation neuronale des moyens de comparaison de la figure 10. 
Figure 1 3 : une representation montrant une fonction d'approximation formee de plusieurs droites. 
Figure 14 : deux representations montrant une determination des limites de validite de deux droites 
consecutives. 

Figure 15 : une representation concernant le raccordement de droites consecutives. 

L'invention concerne Papproximation d'une fonction connue uniquement a travers un certain nombre de 
points par exemple Pi... P 6 (figure 1) dans une representation a *deux dimensions. Chacun de ces points 
est defini par un couple de valeurs (x, y) reliant fa grandeur independante x a la grandeur dependante y. 
Par la suite, il sera question d'une paire (respectivement d'un triplet) de points, ce qui met en oeuvre deux 
couples (respectivement trois couples) de valeurs. On ordonne les points d'apres un ordre croissant des 
valeurs d'abscisses X j( ce qui definit un indice i croissant avec lesdites valeurs. Une convention inverse 
peut etre faite en modifiant en consequence I'expose ci-apres. 

Selon l'invention on effectue une approximation de ('ensemble des couples (Xi, Yi), (X 2 . Y 2 ).„ par une 
droite de regression D ayant pour equation : 

(1) 0: y = p.x + q. 



ou x et y sont des variables courantes. 

Pour cela on considere dans une premiere variante trois couples de valeurs, par exemple (X 3 , Y 3 ), (X4, 
Y*), (X s , Y 5 ), et on determine une droite de regression D par equilibrage des ^reurs absoiues. Une erreur 
est mesuree par la difference apparaissant, pour une abscisse x donnee. entre la valeur y du point et 
I'abscisse y mesuree sur la droite de regression. Equilibrer les erreurs sur trois points, consiste a avoir trois 
erreurs egales en valeur absolue avec un signe d'erreur oppose aux deux autres pour le point ayant une 
abscisse x comprise entre les abscisses x des deux autres points. Puis on examine si pour les points 
restants de I'ensemble. I'erreur qui les separe de la droite reste inferieure ou egale en valeur absolue a 
I'erreur prealablement determined pour les trois points selectionnes. Ceci decoule du fait que Ton 
s'interesse a une erreur pire cas relative a I'ensemble de tous les points a prendre en consideration c'est- 
a-dire la plus grande erreur, en valeur absolue, qui existe entre un des points et la droite de regression. Si 
toutes les erreurs sont effectivement inferieures ou egales. la droite est selectionnee pour representor les 
points sinon on recommence les operations avec trois nouveaux couples de valeurs pour determiner une 
autre droite de regression. 

II peut exister plusieurs droites de regression representant tous les points de I'ensemble. Selon la 
methode on determine la droite de regression optimale qui minimise I'erreur pire cas definie prealablement. 

La figure 1 represente un exemple forme de six points Pi a P 6 disposes selon une representation a 
deux dimensions. A des fins d'explication considerons le resultat final. On observe que la droite de 
regression D de la figure 1 est situee de sorte que les erreurs sont egales en valeur absolue pour les points 
P 3 , P* et P 5 . Pour les points Pi, P 2 . P 6 les erreurs sont inferieures aux precedentes. en valeur absolue. 
Dans un ensemble de points P, a P 6 , la methode va ainsi consister a rechercher les trois points 
particuliers. ici P 3 . Pi, P 5 . qui permettent de determiner la droite de regression optimale qui minimise 
I'erreur pire cas puis a coder cette droite. Dans fe cas de la figure 1 representant un resultat final, si on 
trace deux droites D, et D 2 , paralleles a la droite de regression D, qui passent respectivement par les 
points P 3 , P 5 , d'une part et P* d'autre part, on constate que tous les points de I'ensemble sont a I'interieur 
d'un bandeau limite par les droites Di et D 2 ou sur ces droites. 

La phase de determination de la droite de regression peut donner lieu a plusieurs mises en oeuvre dont 
seules les plus avantageuses seront decrites ci-apres. 
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Premiere vanante de la premiere phase de la methode . 

La figure 3 represente la suite des etapes a mettre en oeuvre pour determiner la droite de regression 
optimale. 

5 Parmi ('ensemble des points Pi... P N (bloc 100), on selectionne (bloc 102) trois points quelconques P,. 
P j( Pk avec i < j < k. Ces trois points servent a determiner la droite de regression qui minimise l'erreur sur 
y pour ces trois points. Cette determination est effectuee de maniere analytique preference I lement par des 
moyens programmes. On determine une droite de reg.ession D pour que trois erreurs (bloc 104) : 
Epo (P if D), E PD (P it D), Ep 0 (P kt D) 

w entre la droite D et chacun des points verifient : 

Epo (Pi. D) = -E P0 (Pj. D) = Ep 0 (P k .D) 
avec 

Epd (Pi. D) = Y; - (p.Xj + q) 

et des relations analogues pour les autres erreurs. 

On determine la droite de regression D a I'aide des coefficients p et q de ('equation (1) tels que : 

20 




25 



XfYj - Xj.Y, + X k .Yj - X,.Y k + X k .Y ± - X, . Y k 
30 2 (X k -X i ) 

L'erreur se rapportant a un triplet (P it P j( P k ) s'ecrit alors : 
35 E T (P i ,P i .P k ) = \£pDiP h D)\ 

Lorsque l'erreur Et a ete ainsi calculee, on examine si les autres points de I'ensemble engendrent des 
erreurs inferieures ou egales. en valeur absolue, a celles des points P i( P j( P k . Pour cela on selectionne un 
point P m (bloc 106) additionnel et on calcule ta valeur absolue de l'erreur E Pm (bloc 108) apparaissant entre 

40 la valeur de la variable y au point P m et la droite D. 

Lorsque cette erreur E Pm est inferieure ou egale en valeur absolue a E T (bloc 110) (repere Y). le point 
P m additionnel est accepte et la methode se poursuit (bloc 112) avec un point additionnel suivant (bloc 
106). Si tous les points additionnels satisfont au critere | E Pm \ S E T , ia droite D est acceptee et ses 
coefficients sont utilises pour coder la droite de regression optimale D op , (bloc 114). 

45 Lorsque cette erreur E Pm est superieure a E T (bloc 110) (repere N), le triplet de points P it P j( P k 
selectionne n'est pas accepte et un autre triplet de points (bloc 116) est choisi (lien 101) dans I'ensemble 
de points (bloc 102). La methode s'acheve a I'obtention de la droite D (bloc 114) satisfaisant ce critere 
meme si tous les triplets de points n'ont pas ete examines. 

II peut apparaitre des situations ou a Tissue de I'etape 116 tous les points possibles ont ete examines 

so et aucun des triplets n'a fourni de solution (bloc 118). Dans ce cas, il est possible de reprendre le 
deroulement de la premiere variante en rempla^ant le test du bloc 1 1 0 par le test suivant : 

\E Pm \S a.E T 

ss ou or est un coefficient legerement superieur a 1. Dans cette meme situation, il est egalement possible de 
faire appel a la deuxieme variante. 

La premiere phase peut comprendre les etapes suivantes : 
A - selection de trois couples de valeurs parmi ladite suite. 
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B - calcul de la fonction lineaire courante de regression D et determination dune erreur de triplet 

Et-\E pd \ s'y rapportant, 
C - selection d'un couple additionnel, 

D - calcul d'une erreur additionnelle E Pm entre le couple additionnel et ladite fonction. 

E - et lorsque | E Pm \ £ E■t{\^Q) pour le couple additionnel. la methode reprend a I'etape C avec un 

couple additionnel suivant, 
F - et lorsque | Ep m \ > E T pour au moins un couple additionnel, 

la methode reprend a I'etape A avec une selection d'un autre groupe de trois couples parmi ladite suite 
de couples, 

G - et lorsque | E Pm \ < E T pour tous tes couples additionnels, la fonction lineaire courante de 

regression est codee et stockee en tant que fonction lineaire d'approximation. 
On peut choisir de scruter I'ensemble des triplets en prenant un ordre croissant ou un ordre decroissant 
ou un ordre aleatoire pour effectuer cet examen. Le triplet qui sera retenu pour determiner la droite de 
regression pourra de ce fait etre detecte a un instant quelconque du deroulement de cette scrutation. II 
s'ensuit que la rapidite d'obtention de la droite de regression depend de I'instant au cours duquel le triplet 
est detecte. Sa mise en oeuvre presente un degre de complexite allant de N a N* ou N est le nombre de 
points initiaux. Sa complexite est done reduite pour un petit nombre de points. Cette variante permet 
d'obtenir une realisation materielle avec une forte parallelisation. Elle est tres peu sensible a une troncature 
des valeurs et fournit un resultat exact. 

Deuxieme variante de la premiere phase de la methode. 

Dans cette deuxieme variante (figure 4), on selectionne successivement des triplets de points, on 
calcule chaque fois une droite de regression, et, par recurrence, on selectionne celle qui delivre I'erreur E P0 
la plus grande e'est-a-dire correspondant a I'erreur pire cas pour I'ensemble de points considered. 
Selon la seconde variante, la premiere phase comprend les etapes suivantes : 
A - selection de trois couples de valeurs parmi ladite suite, 

B - calcul d'une fonction lineaire courante de regression D et determination d'une erreur de triplet 
E T = |Epo| s'y rapportant, 

C - comparaison de I'erreur E T avec une erreur optimale E op initia'isee a une valeur strictement 
negative, 

D - et lorsque E T > E op . mise a jour de I'erreur optimale E op en remplagant Eo P par E T et mise a jour 
des codes d'une fonction lineaire optimale de regression D op en remplagant ceux-ci par les 
codes de la fonction lineaire courante D, 

E - puis retour a I'etape A pour selectionner trois autres couples, 

F - et lorsque tous les triplets de couples de valeurs de la suite ont ete testes, les derniers codes de 
la fonction lineaire optimale D op constituent les codes de la fonction lineaire d'approximation. 

A chaque examen d'un triplet, on compare I'erreur E T dudit triplet avec I'erreur optimale precedemment 
memorised et on met a jour I'erreur optimale E op avec la plus grande valeur de I'erreur E T determined pour 
chaque triplet. On met a jour egalement les parametres de la droite optimale. Avant la mise en oeuvre de 
I'etape A, il taut initialiser la valeur E op a une valeur faible negative, par exemple -1. 

Dans ce cas. la rapidite d'obtention de la droite de regression ne depend pas du mode de scrutation 
des triplets. Sa complexite est de degre N 3 . Une mise en oeuvre materielle peut beneficier de la grande 
regularity de Talgorithme mis en oeuvre. II est peu sensible aux troncatures des donnees et fournit une 
solution exacte. 

Troisieme variante de la premiere phase de la methode. 

Dans cette troisieme variante (figure 5), on selectionne d'abord une paire de points a laquelle on ajoute 
un point supplemental, situe entre ces deux points, afin de former un triplet de points. Pour cela. on 
modifie les etapes A, B et C de la premiere variante, les autres etapes restant les memes. Les etapes 
modifiees sont telles que : 

A1 - modifie I'etape A en operant une selection de deux couples de valeurs appartenant a ladite 
suite, tel qu'il existe au moins un couple additionnel intermediate ayant une grandeur indepen- 
dante (Xt - Xe) comprise entre les grandeurs independantes dudit couple pour constituer au 
moins un triplet de couples, 

A2 - modifie I'etape A premierement en determinant une fonction lineaire annexe qui contient les 
deux couples selectionnes et deuxiemement en determinant des secondes erreurs entre les 
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grandeurs dependantes des couples intermediates possibles et ladite fonction lineaire annexe : 
et lorsque ces secondes erreurs sont toutes de meme signe. selection du couple 
intermediaire fournissant la plus grande seconde erreur, en valeur absolue. pour former un 
triplet de couples de valeurs forme du couple intermediaire et des deux couples 
selectionnes, 

«rt lorsque ces secondes erreurs sont de signes differents. reprise de la methode a I'etape 
Bl - I'etape B est effectuee avec ledit triplet selectionne. 

Ci - modifie I'etape C en selectionnant un couple additionnel dont la grandeur independante n'est 

pas comprise entre la grandeur independante des deux couples selectionnes. 
Lorsque I'erreur | £ Pm | est superieure a I'erreur E T (bloc 110}, la methode reprend a I'etape A1 (lien 
101) avec une nouvelle selection de paire de points (bloc 102a). 

On observe que le deroulement de cette troisieme variante depend de la scrutation des valeurs et done 
des valeurs elles-memes. La complete de la mise en oeuvre de cette variante varie entre N et N 3 ce qui 
lu, confere un certain avantage par rapport aux variantes precedentes. Une forte parallelisation des moyens 
de m.se en oeuvre peut etre operee mais I'implementation des moyens peut presenter un manque de 
regulante ce qu, peut constituer un handicap pour une realisation integree. Cette variante est peu sensible 
aux erreurs d arrondis des valeurs et fournit une solution exacte. 

20 Quatrieme variante de la premiere phase de la methode. 

Elle concerne la determination de la droite de regression a partir des enveloppes 

II est possible de reduire le nombre de triplets a examiner en determinant des enveloppes respective- 

,wlr h rr e e ! ,n , erieUre em ° Urant 165 eXtr§meS dans la presentation bidimensionnelle de 

I ensemble de points. Une enveloppe superieure ou une enveloppe inferieure est definie telle qu'en joignant 

d6 " X P ° imS adiaC6ntS ^ elco ^^s de I'enveloppe, tous les autres points soient situes d'un 

ZZT»n^ t '' enve, °PP e Sl "* rieure ° u * "'enveloppe inferieure. On determine ainsi tous 

les points appartenant a ces dites enveloppes. 

dW^2T inati ° n ^ ' a dr ° ite de ^ 9reSSi0n Va COnsiSter a consid§rer ,es P aires & Points adjacents 
d une des enveloppes auxquels on associe un point intermediaire n'appartenant pas a ladite enveloppe 
pour constituer un triple, et operer comme cela vient d'etre decrit dans le cas des paires de points de to 

e reTopp e e Var,an,e - Si S0 ' U,i0n ° Ptima,e n ' 3 PaS ** **"*»• °" Consid - e ,es ** 

I'envebooT On XT* T Tf PP6 ' ° n S<§ ' eC,ionne une P aire de P° ints Scents appartenant a 
l enveloppe. On determ.ne alors s'H existe un point dispose de telle facon que son abscisse soit 

autre paire de pants de la meme enveloppe. Pour certain* paires. lorsqu'il apparait qu'il existe un ou 

conZnttr poin,s , j rr diaires - on choisit ,e poim ■• Z Ti «2 

con.enant les deux points de la paire pour former un triplet e. pour determiner une droite de regression 
" 51 ce,te droi,e de "*™*n peut etre selectionnee comme droite de regr l St' 

toXeZZe" 0SUVre m§meS 0P§rati ° nS PUe Ce,teS P-'^men, dans le cas de fa 

nrJS» °n? ° n H° d ! fie ,a . ,roisi ° me variante telle <** (figure 6), prealablemen. a I'etape 102a (figure 5) la 
inSZe rl' a ;T C ° mPrend ^ (b ' 0C 100a) de ^terminaton d'une eTeloppe 
S ^ a slction Z T SUPeneUre r4uniSSan ' ' eS POintS ,es P'^ extremes de .'ensemble de 
Zl ou ^TZlirT , P0,mS 3 ''' taPe 1023 ' aite Parmi tes P° ints ad i ace " ts appartenant a 
ooint tnJrml a PP eCt '° n d ° ' adite P3ire de p0ints est effectu * e lorsqu'il existe au moins un 

point intermediate ayant une abscisse situee entre les abscisses des points de .a paire de points S'H 

droite passant par les deux po.nts qu, forment la paire de points. Si une solution n'est pas trouvee avec la 
premiere enveloppe. on poursuit le traitement avec la seconde enveloppe 

fiqureVftdTS T* * SUP * rieure est ef,ectu6e selon ''organigramme de la 

X^ToT^ ^ ,n,erieUre ' ' eS P ° in,S Wf-^ - ''enveloppe inferieure sont 
reperes par la lettre Q. Un point courant Q v est repere" par I'indice v 

coullTeVrepirl T- T* S ° m : °' = P » et Q ' = P " Le «*™< ^ 

courant a v est repere . » - i p0U r P, Le nombre k de points traites est comptabilise. 
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Etape 402 : un premier test est effectue pour determiner si le dernier point P N -i a ete traite pour 
detecter la fin de la determination de I'enveloppe. 

Etape 404 : dans le cas contraire on teste si v > 1 . Lorsque v < 1 on incremente v (v = v + l ) et on 
prend le point P k comme point Q v (etape 407). On incremente rindice k pour traiter le point P suivant 
5 (etape 409) et reprise de la methode a I'etape 402. 

Lorsque v £ 1 , on calcule une droite passant par les points Q v - , et Q v (etape 406) et on determine le 
signe t de I'erreur sur la variable dependante entre le point P k courant et cette droite (Q v -i, Q v ). Ceci a 
pour but de determiner si le point courant est au-dessus ou en dessous de la droite (Q v _ Q v ). 

Lorsque te signe e £ o. il faut supprimer le dernier point Q V( decrementer v tel que v = v - 1 (etape 
io 410) et reprendre la methode a I'etape 404. On peut de cette maniere etre conduit a supprimer certains 
points deja accepted lorsqu'un point suivant necessite de les supprimer. 

Lorsque le signe * est strictement positif, ta methode reprend a I'etape 407 avec un point suivant. 
Get organigramme s'applique a la determination des enveloppes inferieure et superieure en inversant le 
signe de I'erreur a considerer. 
/s Pour faire comprendre les mecanismes ainsi mis en oeuvre considerons. a titre d'exemple. le cas 
simple forme par les points Pi, P 2 , P3. P4 de la figure 2 et determinons I'enveloppe inferieure. Le point Pj 
est le premier point de I'enveloppe d'ou Q Q = Pi . On prend P2 comme deuxieme point Q1 . On calcule la 
droite Lo reliant Pi, P 2 . Le segment Pi, P 2 est accepte comme premier segment de I'enveloppe basse. On 
examine si P 3 est au-dessus de la droite L^, ce qui est le cas sur ta figure 2. Le segment P 2 P 3 est accepte 
20 comme deuxieme segment de I'enveloppe basse. On calcule ensuite la droite Li reliant P 2 et P 3 et on 
examine si P4 est au-dessus de la droite Li . Ce n'est pas le cas sur la figure 2. Ceci entraTne que le point 
P 3 et la droite Li sont d'abord ecartes. De meme, P4 est en dessous de la droite U done le point P2 et la 
droite !_<, sont ensuite aussi ecartes comme ne pouvant pas appartenir a I'enveloppe inferieure et les 
operations se poursuivent avec la determination de la droite L 2 reliant Pi, P 4 , puis I'examen des points 
25 suivants (non represented). 

La complexity de mise en oeuvre de cette variante basee sur des enveloppes est de degre N 2 et 
depend de I'ordre de scrutation des donnees. Cette complexity est moindre que celle des variantes 
precedentes et de ce fait fournit un resultat rapidement. La regularite de I'implementation est moyenne mais 
cette variante est peu sensible aux erreurs d'arrondi des valeurs et delivre une solution exacte. 

Pour certaines applications, i' peut etre souhaitable d'accroltre la precision de la determination de la 
fonction d 'approximation dans certains domaines de la grandeur independante x et d'affecter des coeffi- 
cients de ponderation Wj aux points Pj. C'est par exemple le cas lorsque la fonction d 'approximation est 
faiblement variable avec la grandeur independante x. A certains points peuvent alors etre affectes des 
coefficients de ponderation. Ceux-ci peuvent etre communs a plusieurs points ou etre individuels pour 
chaque point. Ces coefficients de ponderation sont, par la suite, considered strictement positifs. 

Dans les cas ou des coefficients de ponderation existent, on definit une erreur E PD entre un point Pj et 
la droite de regression D tel que : 

E PD (P»D) = W,\Y,-{pX + Q)] 

ou E P0 est une valeur signee. La determination de la droite de regression D pour trois points P jt Pj, P k est 
alors modifiee en ce que les quantites p et q de I'equation (1) deviennent : 

(2) p = NUMPOET et q = NUMQ/DET 

ou les quantites NUMP, NUMQ et DET sont definies par : 

DET = W.W^XpXi) + Wj.WK^X.-Xj) + W^k^X.-XO 
NUMP = Wi.W^Yj-Yj) + Wj.W.^YrYj) + W^W^Y,) 
NUMQ = W i .W j .<X,.Y r X i .Y j ) + W^X^-X^) + W,.W k .(X k .Y r X,.Y k ) 

Par ailleurs. I'erreur E T associee a ce triplet peut etre exprimee et calculee par : 
avec 

ETiPi.P h P k ) = \E PD (P h D\ 
Ep 0 {Pi,D) = -£p 0 (P /( D) = E PD {P k ,D) 

ou 
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(3) P - ^^^\lX i .Y r X 1 .Y i )^(X j .Y k -X k .Y i )*(X k .Y^ 
7 DET 

5 

La premiere et la seconde variante de ta methode de la premiere phase decrites precedemment 
peuvent etre mises en oeuvre en faisant intervenir les coefficients de ponderation ci-dessus. Cette mise en 
oeuvrc peut etre effectuee en programmant un catculateur. 

Dans le cas ou Ton selectionne une paire de points en faisant intervenir des coefficients de ponderation 
io W. la methode presente les amenagements suivants. 

Pour les points P i( P,, P k qui forment le triplet, on definit un facteur F, qui lie les coordonnees 
respectives et les coefficients de ponderation respectifs des points P it P,, P k . Au point central 1 , on associe 
le facteur F, tel que : 




20 Ce facteur F, intervient dans la determination du point intermediaire P ( qui sera a selectionner (i < I < k) 
pour constituer le triplet en tenant compte des poids affectant les points. On choisit un point Pi et on 
calcule la droite de regression Di associee au triplet P j( P t , P k et I'erreur E T i associee au triplet. 

Cette variante de la methode consiste a modifier uniquement I'etape A2 de la troisieme variante (bloc 
102b, figure 5). Cette etape determine I'existence et la valeur d'un point intermediaire servant a former un 
25 triplet. On cherche d'abord a former une droite de regression situee en dessous de Pi et de P k . Pour 
chaque point intermediaire P, f on determine si F, = 1 et E P0 (P it D,) < 0. Si au moins un point veVifie cette 
condition, il n'existe pas de droite de regression situee en dessous de P ; et P kl sinon on determine une 
grandeur G^* qui est la quantite maximale parmi : 

. d'une part les quantites E P0 (P it D p ) pour tous les points intermediates, 
30 . d'autre part les quantites (a) suivantes, uniquement pour les points intermediaires pour lesquels F, < 1 
avec 



35 < 3 > T~^\- EpD (P i'°i> • 

S*il extste au moins un point intermediaire tel que F, > 1. on definit egalement une grandeur G min qui 
est la valeur minimale prise par (a) uniquement pour les points intermediaires tels que F, > 1. On examine 
40 alors s'il existe au moins un point intermediaire tel que : 
Epd (P i( D,) 2 0 

G m ax S a En 

et G min . a Z En (a, coefficient £ 1). S'il existe un tel point, il est choisi comme point intermediaire pour 
former le triplet. 

45 S'il n'existe pas de point intermediaire tel que F, > 1, on teste s'il existe au moins un point 
intermediaire tel que : 
Epd (Pi. Di) 2 0 

G max S or-En (o : coefficient 2 1). 

Si un tel point existe, il est choisi comme point intermediaire pour former le triplet, 
so Si aucun triplet n'a ete constitue, on cherche a former une droite de regression situee au-dessus des 
points Pj et P k . La meme methode est reprise en inversant le signe des erreurs E P0 . 

Si aucun point P, ne peut etre selectionne, on recommence avec une autre paire P i( P k . 
Lorsque I'ensemble de points a traiter est trop important pour etre represente par une seule droite de 
regression, la methode code alors plusieurs droites de regression chacune etant determinee saton la 
55 methode decrite precedemment. 

La figure 13 represente un exemple dans lequel la fonction d'approximation est constitute de plusieurs 
droites de regression. 
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Dans une premiere situation, de par la connaissance des grandeurs a trailer, on peut vouloir imposer 
des limites a chaque droite. de regression suivant la grandeur independante x. Ainsi on peut vouloir 
disposer d'une droite de regression D a entre les valeurs [X a ,XJ de la grandeur, borne x a incluse, borne x d 
exclue. De meme avec D b pour [X b , X c [ et D c pour [X c . X^ . Dans ce cas le probleme revient a determiner 
une droite dans un domaine limite et a appliquer a chaque fois la methode deja decrite. 

Mais il est possible de faire que pour chaque droite. la methode determine des limites suivant la 
grandeur independante x sans qu'eiles soient fixees au depart. 

Le principe de la determination d'une limite optimale entre deux droits de regression adjacentes est 
represents sur la figure 14-A. Soient deux droites de regression D1 et D2 non optimales. La droite D1 est 
determinee a partir de N1 points et la droite D2 est determinee a partir des N2 points restants avec N1 + 
N2 = N, ou N est le nombre total de points. Les droites D1 et D2 torment une approximation avec 
respectivement des erreurs maximales E1 et E2. Supposons que E2 soit superieure a E1 et que Ton 
cherche a reduire E2, il faut alors transferer vers D1 un point appartenant a D2. Lorsque le nombre de 
points diminue. I'erreur resultante est soit inchangee soit reduite. Lorsque le nombre de points augmente, 
I'erreur resultante est soit inchangee soit augmentee. En fonction du nombre n de points appartenant aux 
droites D1 et D2, les variations des erreurs El et E2 peuvent presenter les allures representees sur la figure 
14-B pour I'exemple choisi. II s'ensuit que ('approximation obtenire globalement par D1 et D2 est optimale 
lorsque ces erreurs sont voisines, c'est-a-dire situees dans la zone delimitee par un trait discontinu sur la 
figure 14-B. Ces erreurs ne sont pas necessairement egales compte tenu du caractere discontinu des 
valeurs de la grandeur independante. 

Pour determiner la valeur d'abscisse X lirT1 entre les deux droites : 

- on determine la droite D1 sur un certain nombre de points et on calcule I'erreur E1 maximale. 

- on determine la droite D2 sur les points restants et on calcule I'erreur E2 maximale. 

- on compare E1 et E2 et on transfere un point de la droite qui presente la plus forte erreur vers la 
droite qui presente la plus faible erreur, 

- on determine la valeur limite X Hm lorsqu'il se produit une inversion dans le rapport entre les deux 
erreurs. 

Cette methode s'applique par recurrence a plusieurs droites dSlimitees entre elles par plusieurs limites 
d'abscisses. 

Les droites sont determinees a partir d'une suite discrete et limitee de grandeurs de mesure. 
Neanmoins, pour I'exploitation des droites de regression, il est necessaire de definir leur domaine 
d'existence qui s'etend sur un continum de valeurs situees entre deux limites d'abscisses. Or ia determina- 
tion des droites fournit une suite de droites qui ne se raccordent pas necessairement par leurs extrSmites. II 
peut apparaitre utile pour certaines applications d'eviter qu'un ecart apparaisse sur la valeur de la grandeur 
dependante y pour des valeurs voisines (X lim - <) et (X lim + «) de la grandeur independante ou « est une 
tres faible valeur. tl est possible de choisir que I'abscisse limite X| im appartienne exclusivement a I'une ou a 
Tautre droite. II est aussi possible de faire suivre la determination des suites de droite d'une procedure de 
raccordement de droites. 

Ceci est represents sur la figure 15. Une solution peutconsister a valider la droite D1 jusqu'a I'abscisse 
du premier point appartenant a la droite D2 et a recalculer ia droite D2 a partir de la nouvelle valeur de 
grandeur dependante y ainsi determinee. En conservant le dernier point appartenant a D2, on determine 
ainsi une nouvelle droite D'2 representee en pointilte sur la figure 15. Une procedure analogue peut etre 
appliquee pour substituer la droite D'3 a la droite D3. On obtient ainsi un ensemble forme de plusieurs 
droites de regression forme, dans I'exemple decrit, par les droites D1, D*2, D'3. Cet ensemble constitue une 
approximation des grandeurs de mesure en reduisant une erreur maximale entre les grandeurs de mesure 
et ('ensemble de droites. Une variante plus adaptee a cette reduction consiste a faire que ce soient les 
droites correspondant aux plus fortes erreurs qui imposent leurs points limites comme points limites aux 
droites correspondant a des erreurs plus faibles. 

Dispositif cour une mise en oeuvre de la premiere phase 

Un schema d'un dispositif pour generer une fonction d'approximation selon I'invention est represents 
sur la figure 8. Le dispositif 5 regoit des couples (X;. Yj) associant la grandeur dependante Y 4 a la grandeur 
independante X f . Les couples (X j( Yj) entrent dans les premiers moyens 10 pour determiner et coder la 
fonction lineaire de regression formant une approximation des couples de mesure. Les codes specifiques 
ainsi determines sont transmis (connexion 9) aux seconds moyens 17 qui determined des seconds couples 
(X A . Y* a ) d'apres les codes specifiques a la fonction d'approximation. 
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On distingue d'une part la seconde variante qui determine toutes les droites de regression passant par 
chaque combinaison realisable de triplets avec la suite de premiers couples, et les autres variantes pour 
lesquelles a chaque droite de regression determinee (associee a son erreur E T ). on examine si les points 
additionnels restants fournissent bien des erreurs additionnelles inferieures ou egales a I'erreur de triplet E T . 

La figure 9 represente un dispositif adapte a la seconde variante. II comprend : 

- une memoire 12c MEM qui stocke notamment tous les points appartenant a I'ensemble de points a 
traiter. Ces points sont representes par leurs coordonnees (x. y). et eventuellement leurs poids W ou 
leurs inverses 1/W, 

- une unite de calcul 13c COMPUT qui calcule. pour chaque triplet selectionne. la droite de regression 
adaptee a chaque triplet, c'est-a-dire les codes p. q de la droite et I'erreur E T associee au triplet 

De plus, un controleur 11c CONTR permet d'effectuer la gestion des operations et d'adresser des 
nouveaux triplets en effectuant les lectures/ecritures de la memoire 12c et le chargement de I'unite de 
calcul 13c par de nouveaux triplets. La selection de la droite de regression qui est a conserves c'est-a-dire 
pour cette variante celle delivrant la plus grande erreur de triplet, est effectuee par l'unite de calcul 13c. 

La figure 10 correspond au cas des autres variantes pour lesquelles pour chaque droite courante de 
regression on examine si les autres points additionnels de la suite dSlivrent une erreur inferieure a I'erreur 
de triplet. 

Pour cela. les premiers moyens 10 comprennent une unite de calcul 13c et une unite de comparaison 
14c formant des moyens de calcul 19c. 

L'unite de calcul 13c. transmet les codes p. q de la droite de regression du triplet courant a I'unite de 
comparison 14c COMPAR qui determine si les points additionnels de ('ensemble de points generent une 
erreur plus fable avec cette droite de regression que celle generee par les points du triplet courant Pour 
cela, l'unite de comparaison 14c effectue le test : 

I E Pm \> E T 

Si la droite de regression courante n'est pas acceptee (test positif). un autre triplet est selectionne et une 
procedure analogue est a nouveau effectuee. Si le test est negatif pour tous les points additionnels la 
droite de regression est acceptee et ses parametres sont charges par I'unite de calcul 13c dans la memoire 

Pour les variantes qui operent a partir d'une paire de points. I'unite de calcul 13c est programmee pour 
determiner les points intermediates et pour former les triplets issus des paires de points, ceux-ci etant eux- 
memes eventuellement issus d'enveloppes de points. Dans ce dernier cas. I'unite de calcul 13c est 
egalement programmee pour determiner les enveloppes. L'unite de comparaison 14c effectue ensuite les 
comparisons des erreurs se rapportant aux differents triplets. 

Pour effectuer la gestion des operations et adresser des nouveaux triplets (lorsque des triplets 
examines nont pas donne satisfaction), le dispositif comprend. comme precedemment. un controleur 11c 
QUI r© . 

- les lectures/ecritures de la memoire 12c. 

- le chargement de I'unite de calcul 13c par de nouveaux triplets. 

- I'appel vers I'unite de comparaison 14c de tous les points additionnels a examiner successivement 

II est possible de mettre en oeuvre la premiere phase de la methode a I'aide de moyens de 
dTSSl h 9UrC H- 2) ayam StfUCtUre neur0nale ' Consid e™s «e cas le plus complet ou il existe 

?ieur doTner V'T'^l'T^ * "** * L ° W ^ COeffiCientS n ' 8XiS,ent »** « su «« 

TJZ IT To V f, Uf Un ' te d3nS ' eS ex P ,ications 0" v °"« suivre. A titre d'exemple. on peut stocker 

J un te de ^r? f h 9 ' 10> 165 i0VerSeS dSS COeff ' Cien,S Wi af,ect « '"dividuellement aux points P,. 
L unite de calcul 13c determine selon les equations (2) et (3). pour un triplet donne, les codes : -p. q. E T et 

^ La condition a tester dans Tunite 14c pour un point additionnel P m a tester, de parametres X m . Y m . W m , 
Y m - pX m - (E T 'W m ) - q > 0 



ou 



Ym - p.X m + (E r W m ) - q < 0. 



EP 0 624 847 A1 



Ces tests sont aisement realisables par une unite 14c ayant une organisation neuronale. En effet. ce 
test est constitue de fonctions lineaires, de fonctions de seuillage et de fonctions logiques aisement 
realisees dans un reseau de neurones. 

Une telle unite neuronale de comparaison est representee sur la figure 12. Elte comprend trois 
5 neurones N1, N2, N3. Les neurones N1 et N2 regoivent des donnees E1. E2, E3. E4. Le neurone N3 regoit 
les sorties des neurones N1 et N2. Chacune des entrees de ces neurones est affectee d'un coefficient 
synaptique Ct selon la technique connue mise en oeuvre dans les reseaux de neurones. Cette technique 
est par exempfe decrite dans R.P. LIFPMANN, "An introduction to computing with neural nets" IEEE ASSP 
Magazine. April 1987, pp. 4 to 22. 
w Pour effectuer les tests precedents, on programme les neurones N1 et N2 selon le tableau I : 

TABLEAU I 
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Neurone 


E, 


C, 


E 2 


C 2 


E 3 


c 3 


E* 


C* 


1 




1 




" P 




-E T 


1 


-q 


2 


Y m 


1 




-P 


1/w m 


Et 


1 


-q 



20 Tous les coefficients synaptiques du neurone N3 sont egaux a 1 . 

Les neurones N1 et N2 calculent chacun un potentiel neuronal s tel que : s = C,.E, 
A ce potentiel s, il faut appliquer une fonction d'activation A. Pour les neurones N1 et N3 t la fonction 
d'activation Ai est : 

25 At (s) = i si s > 0 
Ai (s) = 0 si s S 0. 

Pour le neurone N2, la fonction d'activation A 2 est : 

30 A 2 (s) = 1 si s < 0 

A 2 (s) = 0 si s £ 0. 

La sortie du neurone N3 vaut 1 si la condition a tester est satisfaite et vaut 0 dans le cas inverse. 

35 On observe sur le tableau I que les donnees -p, -q t -E T , + E T sont des codes qui apparaissent comme 
coefficients synaptiques dans le schema de la figure 12. Les donnees Y m , X m , 1/W m sont des donnees qui 
apparaissent sur les entrees du meme schema. Lorsque tous les points additionnels ont ete testes et 
qu'une droite de regression a ete selectionnee les codes p. q. W charges dans la memoire 12c sont alors 
utilises pour mettre en oeuvre la seconde phase de la methode (figure 11). 

40 L'interet presente par la realisation neuronale decrite est de pouvoir paralleliser les differentes 
operations a mettre en oeuvre selon les variantes deja decrites. Le fonctionnement d'une telle realisation 
neuronale est alors tres rapide. 

Lorsque les droites de regression ont ete determinees, la seconde phase consistant a calculer des 
seconds couples (X A , Y* A ) de valeurs des grandeurs est alors mise en oeuvre dans des moyens de 

45 decodage 17 (figure 11). Les codes des droites sont charges dans une memoire 12a qui au cours de la 
seconde phase est adressee, par un controleur 11 a, pour fournir les codes des droites de regression 
adressees. La memoire 12a, organisee en lignes par exemple. contient pour chaque droite de regression 
les parametres p, q et x L ou x L est la limite superieure d'abscisse pour laquelle chaque droite de regression 
est definie. 

so La memoire 12a contient ainsi un tableau de parametres 

X L1' Pq- q o 



x Lm* p m-l' q m-l 

correspondant aux m droites de regression stockees. 
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Les moyens de decodage 1 7 comprennent : 
!e controleur 1 l a . 
la memoire 12 a , 

et une unite de decodage 13 3 . Celle-ci recoil une valeur de requete X A de la grandeur independante 
pour laquelle on desire obtenir un resultat Y' A issu de la fonction d'approximation. Pour ceta le 
controleur 11 a presente successivement certaines lignes de codes de la memoire 12 a et, pour chaque 
ligne de codes. I'unite de decodage 13 3 teste le code x L pour determiner si la valeur X A de la 
grandeur d'entree peut etre traitee par la droite de regression codee par ladtte ligne. Si cela n'est pas 
le cas. I'unite 13 3 requiert du controleur H a une autre ligne de codes. Lorsque cette condition est 
verifiee. I'unite I3 a effectue le calcul : 

Y'a = pn-X A + q n 

si la ligne d'ordre n est la ligne adequate. 
Comme cela vient d'etre decrit, les premiers moyens 10 peuvent constituer un dispositif de codage et 
les seconds moyens 17 peuvent constituer un dispositif de decodage. les deux dispositi f s etant separes. lis 
peuvent en effet etre distants Tun de ('autre, les operations de codage etant separees des operations de 
decodage. Mais les deux moyens 10 et 17 peuvent aussi former un dispositif de codage/decodage unique. 
Dans ce cas certains organes peuvent etre affectes successivement a plusieurs operations. Notamment : 
20 - il peut s'agir d'une memoire 12 unique qui peut stocker les parametres des couples de valeurs et les 
codes des droites de regression determinees. Les unites 12a et 12c forment alors la memoire 12. 

- il peut s'agir d'un controleur 11 unique qui regroupe les fonctions des controleurs 1 1c et 11a. 

- I'unite de calcul 13c et ('unite de decodage 13a peuvent egalement former une seule unite 13. 

II est egalement possible de realiser les moyens de decodage 17 (figure 11) selon une organisation a 

25 structure parallele. Une solution consiste alors a comparer, en parallele, la grandeur de requete X A a tous 
les codes x L . On obtient un ensemble de signaux T L representant les resultats de ces tests de comparison. 
On realise ensuite. en parallele, pour toutes les valeurs de L, un test qui combine deux a deux les signaux 
T L et T L+1 et qui indique si la grandeur de requete X A est dans rintervalle d'indice L. Ce test est verifie pour 
une unique valeur de L, notee n. On accede alors aux codes p n et q n correspondants, et on en deduit la 

30 grandeur Y' A de la meme maniere que precedemment. 

Les premiers moyens 10 (de codage) associes aux seconds moyens 17 (de decodage) peuvent etr 
utilises pour determiner la valeur d'une fonction d'approximation par au moins une droite de regression. 
Ces determinations peuvent etre effectuees pour des valeurs quelconques de la grandeur independante (a 
Tinterieur de bornes predefines definissant Tetendue de Taction de chaque droite de regression). Cette 

35 methode evite de stocker inuttlement des tables de valeurs, pour lesquelles toutes les valeurs ne seront pas 
exploiters. Selon la methode on ne determine que les seules valeurs necessaires a ('application. La 
methode selon I'invention presente Ttnteret de ne calculer que les valeurs necessaires. Le dispositif de 
Tinvention peut etre utilise en combinaison avec un processeur neuronal pour calculer une approximation 
d'une fonction non lineaire par exemple une fonction sigmoYde. Un processeur neuronal calcule des 

40 potentiels neuronaux issus de sommes de produits de coefficients synaptiques par des etats de neurones. 
La technique est developpee dans le document de R.P. LIPPMANN deja cite. Chaque potentiel neuronal 
doit etre soumis a Taction d'une fonction non lineaire. Pour cela, chaque potentiel neuronal constitue la 
grandeur independante decrite precedemment et le calcul de la fonction th(x) par exemple est effectue pour 
les seules valeurs de requetes utiles. La methode fournit des valeurs approximees, mais cette approxima- 

45 tion peut etre obtenue avec une precision qui peut etre accrue selon le nombre de droites de regression 
calculees et eventuellement les coefficients de ponderation W. La precision reste neanmoins tributaire de la 
precision des coordonnees des couples de depart. 

Une telle methode est particulierement interessante pour les calculs de fonctions connues (telles que 
des fonctions mathematiques) ou de fonctions inconnues explicitement. par exemple une fonction represen- 

50 tee par des points de mesures. que Ton desire simplifier par une fonction lineaire de regression. 

Dans sa version neuronale, Tinvention est interessante dans des applications neuronals car elle 
apporte non seulement un traitement homogene mais aussi une grande compacite de ('architecture 
materielle necessaire. 

II est possible de combiner les vartantes de la methode, dans leurs mises en oeuvre materielles pour 
55 combiner leurs avantages. Ainsi par exemple, les combinaisons suivantes sont avantageuses : 
Dans le cas ou les poids existent, deux combinaisons apparaissent interessantes : 
Dans la premiere combinaison, on utilise tout d'abord la premiere variante basee sur des triplets de 
points. Si une solution n f est pas trouvee. on enchaine avec la deuxieme variante basee aussi sur des 
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triplets de points. Ceci permet d'obtenir les avantages suivants : 
. Quelles que soient les donnees. on obtient une solution. 

Dans le cas ou on exige qu'on obtienne toujours une solution exacte. cette combinaison permet 
d'obtenir la rapidite maximale. En effet. dans le cas general la solution est fournie directement par 
la premiere variante appelee, qui est celle fournissant le plus raptdement une solution exacte dans 
le cas ou les poids existent. La deuxieme variante est beaucoup plus lente mais elle permet de 
garantir que la combinaison fournit toujours une solution. 
. Dans la deuxieme combinaison, on utilise tout d'abord ta variante basee sur des paires de points. Si 
elle ne trouve pas une solution (qui peut etre approchee), on utilise la deuxieme variante basee sur 
des triplets de points. Ceci permet d'obtenir les avantages suivants : 

Quelles que soient les donnees, on obtient une solution (eventuellement approchee). 
Cette combinaison est en moyenne plus rapide que la precedente. 
En contrepartie. elle ne fournit generalement qu'une solution approchee. 
Dans le cas ou les poids n'existent pas, une combinaison interessante consiste a utiliser tout d'abord la 
variante basee sur I'enveloppe des points. Si elle ne trouve pas de solution, on utilise la deuxieme variante 
basee sur des triplets de points. Ceci permet d'obtenir les avantages suivants : 
Quelles que soient les donnees, on obtient une solution. 
Cette combinaison permet d'obtenir la meilleure rapidite moyenne. 
Lorsque la fonction lineaire de regression a ete determined selon la methode decrite, if est possible de 
determiner non seulement des grandeurs Y' A ressortissant a ladite fonction mais egalement des grandeurs 
ressortissant a la derivee de ladite fonction. En effet, pour une valeur de requete X A de la grandeur 
independante, it suffit d'extraire des moyens de memorisation la valeur p ressortissant a (a valeur X A pour 
connaTtre la valeur de cette derivee correspondante. 

Revendic'ations 

1. Dispositif (5) pour generer une fonction d'approximation fondee sur des premiers couples ((Xi, Yt) a 
(Xe, Y 6 )) de valeurs associant une grandeur dependante (Yi a Y 6 ) a une grandeur independante (Xi a 
X6). et pour determiner des seconds couples (X A , Y A ) de valeurs desdites grandeurs d'apres ladite 
fonction d'approximation caracterise en ce que le dispositif comprend : 

- des premiers moyens (10) : 

- pour determiner iterativemant au moins une fonction lineaire courante de regression en 
rendant egales, en valeur absolue. des premieres erreurs, de signes altern^s, mesurees entre 
des premieres valeurs (Y 3 , Y 4 , Y 5 ) de la grandeur dependante pour trois couples (X 3 , Y 3 ) (X4, 
Y4) (X 5 , Y 5 ) d'une suite desdits premiers couples, et respectivement des secondes valeurs 
(Y'3, Y' 4 , Y' 5 ) de la grandeur dependante determinees, d'apres ladite fonction lineaire courante. 
pour les memes valeurs (X3.X4.X5) de la grandeur independante. 

- pour selectionner celle des fonctions lineaires courantes qui delivre I'approximation de tous les 
couples de ladite suite avec des erreurs minimales, 

- et pour coder, a ('aide de codes specifiques (p, q), la fonction lineaire de regression 
selectionnee, 

- et des seconds moyens (17) pour determiner lesdits seconds couples <X A , Y' A ) a I'aide desdits 
codes specifiques. 

2. Dispositif selon la revendication 1 caracterise en ce que les premiers moyens (10) comprennent : 

- des moyens (12c) pour memoriser les grandeurs des premiers couples et lesdits codes specifi- 
ques, 

- des moyens (13c) (19c) pour calculer successivement les codes specifiques (p, q) de fonctions 
lineaires courantes de regression et leurs premieres erreurs. egalisees en valeur absolue. 
respectives (E T ) et pour selectionner au moins une fonction lineaire courante de regression, 

- et des moyens (11c) pour controler les moyens (12c) pour memoriser et les moyens (13c) pour 
calculer. 

3. Dispositif (5) selon les revendications 1 ou 2 caracterise en ce que pour selectionner ladite fonction 
lineaire, les premiers moyens comprennent des moyens (14c) : 

- pour determiner pareillement des erreurs additionnelles pour les autres couples (X lt Yi) (X 2 . Y 2 ) 
(X&, Y 6 ) de valeurs de ladite suite, 

- pour comparer les erreurs additionnelles avec lesdites premieres erreurs egalisees (Et), 
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- et pour selectionner la fonction lineaire courante qui fournit des erreurs additionnelles inferieures 
ou egales en valeur absoiue auxdites premieres erreurs egalisees (E T ). 

Dispositif selon la revendication 2 caracterise en ce que, pour selectionner la fonction lineaire de 
regression, les moyens pour calculer (13c) selectionnent celle des fonctions lineaires courantes de 
regression qui delivre des premieres erreurs egalisees (E T ) maximales. 

Dispositif selon une des revendications 1 a 4 caracterise en ce que les seconds moyens (17) 
comprennent : 

- des moyens (12a) pour memoriser lesdits codes specifiques d'au moins une fonction lineaire de 
regression, 

- des moyens (13a) pour selectionner, sur requete d'une valeur de grandeur independante (X A ), 
une fonction lineaire de regression ressortissant a ladite valeur, et pour calculer, a partir des 
codes specifiques recus des moyens pour memoriser, une valeur de grandeur dependante (Y' A ) 
associee a la valeur de grandeur independante (X A ), 

• et des moyens (11a) pour adresser aux moyens pour calculer, et cjr leur demande, les codes 
specifiques des fonctions lineaires de regression. 

Dans un dispositif pour generer une fonction d'approximation fondee sur des premiers couples ((Xi , Yi ) 
* (X&, Y6)) de valeurs associant une grandeur dependante (Yi - Y&) a une grandeur independante (Xi - 
Xs ), methode pour generer ladite fonction d'approximation a partir des premiers couples (X : , Yi ) (Xg . 
Y & ) de valeurs, et pour determiner des seconds couples (X A , Y' A ) de valeurs desdites grandeurs 
d'apres ladite fonction d'approximation, caracterisee en ce que, la methode comprend : 
une premiere phase : 

• pour determiner iterativement au moins une fonction lineaire courante de regression en 
rendant egales, en valeur absoiue, des premieres erreurs E P0 , de signes alternes, mesurees 
entre des premieres valeurs (Y 3 , Y 4 , Y 5 ) de la grandeur dependante pour trois couples <X 3 , 
Y 3 ), (X4, Y4) (Xs, Y5) d'une suite desdits premiers couples, et respectivement des secondes 
valeurs (Y' 3l Y' 4 , Y' 5 ) de la grandeur dependante determines d'apres ladite fonction lineaire, 
pour les memes valeurs (X3. X*. Xs) de ta grandeur independante, 

* pour selectionner celle des fonctions lineaires courantes qui delivre ('approximation de tous les 
couples de ladite suite avec des erreurs minimales, 

- et pour coder la fonction lineaire de regression selectionnee a I'aide de codes specifiques, 

- et une seconde phase pour determiner lesdits seconds couples (X A , Y' A ) a I'aide desdits codes 
specifiques. 

Methode selon la revendication 6 caracterisee en ce que la premiere phase comprend les etapes 
suivantes : 

A - selection (102) de trois couples (P3. P4. P5) de valeurs parmi ladite suite, 

B - calcul (104) de la fonction lineaire courante de regression D et determination d'une erreur de 

triplet E T =\E P d[ s'y rapportant. 
C - selection (106) d'un couple additionnel (P1), (P2), (P6), 

D - calcul (108) d'une erreur additionnelle E Pm entre !e couple additionnel et ladite fonction, 

E - et lorsque | E Pm | £ £7(1 10) pour le couple additionnel, la methode reprend a I'etape C avec 

un couple additionnel suivant, 
F - et lorsque | E Pm \ > E T pour au moins un couple additionnel (110), la methode reprend a 

I'etape A avec une selection d'un autre groupe de trois couples parmi ladite suite de 

couples, 

G - et lorsque | E Pm | £ E T pour tous les couples additionnels (112), la fonction lineaire courante 
de regression est codee et stockee (1 14) en tant que fonction lineaire d'approximation. 

Methode selon la revendication 6 caracterisee en ce que fa premiere phase comprend les etapes 
suivantes : 

A - selection (102) de trois couples (P3, P4, P5) de valeurs parmi ladite suite. 

B - calcul (104) d'une fonction lineaire courante de regression D et determination d'une erreur de 

triplet £Y = |£pd| s'y rapportant, 
C- comparaison (210) de I'erreur E T avec une erreur optimale E op initialised a une valeur 

strictement negative. 
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D - et lorsque E T > E op . mise a jour (1 14) de I'erreur optimale E op en remplacant E op par E T et 

mise a jour des codes d'une fonction lineaire optimale de regression D op en remplagant 

ceux-ci par les codes de la fonction lineaire courante (D). 
E - puis retour (1 16) a I'etape A pour selectionner trots autres couples (102), 
F - et torsque tous les triplets de couples de valeurs de la suite ont ete testes (116). les derniers 

codes de la fonction lineaire optimale D op constituent les codes de la fonction lineaire 

d'approximation (114). 

9. Methode selon la revendication 8 caracterisee en ce que la premiere phase comprend les etapes 
modifiees suivantes : 

A1 - modifie I'etape A en operant une selection (102a) de deux couples de valeurs appartenant a 
ladite suite, tel qu'il existe au moins un couple additionnel intermediate ayant une grandeur 
independante (Xi - Xe) comprise entre les grandeurs 
independantes dudit couple pour constituer au moins un triplet de couples de valeurs. 

A2 • modifie I'etape A premierement en determinant (102b) une fonction lineaire annexe qui 
contient les deux couples selectionnes et deuxiemement en determinant des secondes 
erreurs entre les grandeurs dependantes des couples intermediates possibles et ladite 
fonction lineaire annexe : 

. et lorsque ces secondes erreurs sont toutes de meme signe. selection du couple 
intermediaire fournissant la plus grande seconde erreur, en valeur absotue, pour 
former un triplet de couples de valeurs forme du couple intermediaire et des deux 
couples selectionnes, 

. et lorsque ces secondes erreurs sont de signes differents, reprise de la methode a 
I'etape A1 , 

B1 - I'etape B (104) est effectuee avec ledit triplet selectionne, 

C1 - modifie I'etape C (106) en selectionnant (106a) un couple additionnel dont la grandeur 
independante n'est pas comprise entre la grandeur independante des deux couples selec- 
tionnes. 

10. Methode selon la revendication 9 modifiee en ce que. prealablement a I'etape A1, lesdits couples de 
valeurs etant represents par des points (P1 - P6) dans un espace a deux dimensions, la premiere 
phase de la methode comprend une etape (100a) de determination d'une enveloppe inferieure et/ou 
d'une enveloppe superteure reunissant les points les plus extremes de I'ensemble de points, la 
selection des deux couples de valeurs a I'etape A1 etant faite parmi les points adjacents appartenant a 
Tune ou I'autre enveloppe. 

11. Methode selon une des revendications 6 a 10 caracterisee en ce qu'a chaque valeur, prise dans la 
suite de valeurs de la grandeur independante, est associee un coefficient specifique de ponderation qui 
pondere lesdites premieres erreurs. 

12. Methode selon une des revendications 6 a 11 caracterisee en ce qu'elte genere une fonction 
d'approximation comprenant une suite de fonctions lineaires de regression raccordees entre elles. 
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